Using the de Haas-van Alphen effect to map out the closed 
three-dimensional Fermi surface of natural graphite 
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The Fermi surface of graphite has been mapped out using de Haas van Alphen (dHvA) measure- 
ments at low temperature with in-situ rotation. For tilt angles 6 > 60° between the magnetic field 
and the c-axis, the majority electron and hole dHvA periods no longer follow the cos(S) behavior 
demonstrating that graphite has a 3 dimensional closed Fermi surface. The Fermi surface of graphite 
is accurately described by highly elongated ellipsoids. A comparison with the calculated Fermi sur- 
face suggests that the SWM trigonal warping parameter 73 is significantly larger than previously 
thought. 



Graphite consists of Bernal stacked graphene layers 
with a weak inter layer coupling which leads to an in- 
plane dispersion which depends on the momentum in the 
direction perpendicular to the layers, k^. Graphite is 
a semi metal with the carriers occupying a small region 
along the H—K—H edge of the hexagonal Brillouin zone. 
The Slonczewski, Weiss, and McClure (SWM) Hamilto- 
nian with its seven tight binding parameters 70, ...,75, A, 
is based on group theoretical considerations and provides 
a remarkably accurate description of the band structure 
of graphite [3 H]. In a magnetic field, when trigonal 
warping is included (73 ^ 0) levels with orbital quantum 
number n couple to levels with orbital quantum number 
n 4- 3 and the Hamiltonian has infinite order. However, 
the infinite matrix can be truncated as the eigen values 
converge rapidly [3 . The validity of the SMW-model, 
has been extensively verified using many different ex- 
perimental techniques e.g. Shubnikov-de Haas (SdH), 
de Haas- van Alphen (dHvA), thermopower, magneto- 
transmission, and magneto-reflectance measurements [H- 
[T5] . However, recently claims [T71 [TH] for the observa- 
tion, in electrical transport measurements, of massless 
two-dimensional (2D) charge carriers with a Dirac-like 
energy spectrum have caused much controversy |19H21j . 

The Fermi surface of graphite has electron and hole 
majority carrier pockets with maximal extremal cross 
sections at ^2 = (electrons) kz ~ 0.35 (holes). For 
both types of charge carriers the in-plane dispersion 
is parabolic (massive fermions). Only at the H point 
{kz = 0.5) the in-plane dispersion is linear, similar to that 
of charge carriers in graphene (massless Dirac fermions). 
At the H-point, there two possible extremal orbits. A 
minimal (neck) orbit of the majority hole carriers, which 
gives rise to minority carrier effects (a-surface) and a 
maximal extremal orbit of the small ellipsoidal minority 
hole pocket (^ surface) which results from the intersec- 
tion of the two majority hole ellipsoids. The existence 
of the three minority carrier pockets at the K-point, the 
so called outrigger pieces, was proposed by Nozieres [22] . 
however their existence is considered to be unlikely due 
to the rather large value of 73, the SWM trigonal warping 
parameter, required. 



In this Letter, we present a complete map of the Fermi 
surface of natural graphite obtained from dHvA mea- 
surements at low temperature (T « 0.4 K) with in-situ 
rotation. For tilt angles 9 < 60°, the dHvA periods of 
both the electrons and holes follow a cos{d) dependence. 
While such a quasi-2D behavior is well established in the 
literature, previous dHvA measurements [S], were unable 
to distinguish between a highly elongated 3D ellipsoid 
and a cylindrical 2D Fermi surface. Our results at larger 
tilt angles demonstrate unequivocally that graphite has a 
3D closed Fermi surface which is accurately described by 
highly elongated ellipsoids provided the spin splitting is 
included. A comparison of our data with the full SWM 
calculations allows us to refine the SWM tight binding 
parameters, notably 73 is found to be significantly larger 
than previously thought. 

For the dHvA measurements we used a mm-size piece 
of natural graphite, which was mounted on a CuBe can- 
tilever which forms the mobile plate of a capacitive torque 
meter. The capacitive torque signal was measured with a 
lock-in amplifier using conventional phase sensitive detec- 
tion at 5.3 kHz. The measurements were performed us- 
ing a 16 T superconducting magnet and a dilution fridge, 
equipped with an in situ rotation stage. Fig. [I] (a) shows 
the torque r as a function of the total magnetic field 
from B = 0- 0.21 T for a tilt angle 6* = 16°. The torque 
shows the expected dependence, t{B) oc —B^ (broken 
line), since r = MBsm{9) and the magnetization M de- 
pends linearly on the magnetic field. Superimposed on 
the large monotonic background, small quantum oscilla- 
tions are clearly visible, which reflect the oscillatory mag- 
netization of the system as Landau levels pass through 
the Fermi energy. The dHvA oscillations can be bet- 
ter observed in the oscillatory torque {tosc)- Here the 
monotonic background has been removed by subtract- 
ing a smoothed (moving window average) data curve. 
The torque t oc sin{9) so that we cannot directly ac- 
cess the quantum oscillations in perpendicular field. In 
order to compare with our previous magnetotransport 
measurements we use the low angle 9 = 16° data writing 
B_L ^Bcos{9). Fig.[l];a) shows Tosc{B). 

The phase and the frequency of the dHvA oscillations. 
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FIG. 1. (color online) (a) The torque r and oscillatory torque 
Tosc versus total magnetic field for 9 — 16°. (b) Fourier 
transform of r(l/Bcos(0)). (c-d) The phase shift function 
5R [exp (193) /(Bx)] as a function of the phase and frequency, 
(e) Oscillatory torque Tosc{B) for B \\ ah {9 = 90°). 



were extracted from a Fourier analysis. The fundamental 
frequencies = 6.32 ± 0.1 T and = 4.62 ± 0.1 T 
for electrons and holes respectively, are obtained from 
the amplitude of the Fourier transform of Tq^c 
(see Fig. [l] (b)). In Fig. [IJc) and (d) we plot the 
phase shift function K{ip,B±) = ^[exp {i(p) f{B±)] — 
cos{ip—(fo)F{B±) as a function of the perpendicular mag- 
netic field and the phase. The fundamental frequency 
and the phase ipo can be found from the maxima in the 
(fi — B plane. The phase values in units of 27r obtained 
are (/sg = -(0.48 ± 0.1) and (p^ = -(0.56 ± 0.1). The 
phase (po = 7 — (5 with 7=1/2 for massive Fermions 
or 7 = for massless Dirac fermions. For a 3D Fermi 
surface the curvature along gives S = ±1/8 for min- 
imum/maximum extremal cross sections. In contrast a 
cylindrical 2D Fermi surface gives S — 0. We can there- 
fore conclude, in agreement with recent dHvA measure- 
ments on highly oriented pyrolytic graphite (HOPG) [l5] . 
that both the electrons and holes are massive Fermions 
with a parabolic energy spectrum (i.e. 7 = 1/2). In 
Fig. [T|e) we show the oscillatory torque Tosc{B) in the 
B\\ab configuration [9 — 90°). The B \\ ab configuration 
can be found very precisely [59 < 0.1°), since the magne- 
tization background changes sign at 6* = 90°. Well pro- 
nounced quantum oscillations are observed demonstrat- 
ing unequivocally that the Fermi surface of graphite is 
3D and closed. 

In order to map out the Fermi surface, we have per- 
formed systematic angle dependent measurements. In 
Fig. [2|^a) we plot the amplitude of the Fourier transform 
of Tosc{^/B) as a function of the period of the oscilla- 
tions) and the tilt angle 9. The hole and electron fea- 
tures, together with the hole harmonic, can clearly be 
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FIG. 2. (color online) (a-b) Color plot : Amplitude of the 
Fourier transform of Tosc{^/ B) as a function of the period of 
the oscillations and the tilt angle (9). The calculated dHvA 
periods for ellipsoidal electron and hole Fermi surfaces are 
also shown for hole (solid lines), hole harmonic (thin solid 
lines), and electrons (dot-dash). The hole neck orbits with 
and without breakdown are shown as dotted lines. Woollam's 
minority carrier data[23] is shown as symbols. The dashed 
line is the prediction for an ellipsoid, (c) Woollam's minority 
carrier data plotted over an extended range. The prediction 
for an ellipsoid (dashed line) and for a neck orbit (dotted line) 
are shown for comparison. 



distinguished. For angles 9 < 60°, the dHvA period for 
both electrons and holes follow the well documented [51 |S] 
cos(6') dependence. Such a behavior is characteristic of 
either a 2D material or very anisotropic material with an 
almost perfectly cylindrical Fermi surface. We can dis- 
tinguish between these two scenarios at higher tilt angle. 
The non cylindrical nature of the Fermi surface is clearly 
revealed for 9 > 60° where deviations from the cos{9) 
behavior, are observed. Namely for 9 > 75° (see inset of 
Fig. [2]), the slope of the dHvA periods for both the elec- 
trons and holes features changing dramatically reaching 
almost zero close to 9 = 90° (see Fig.[2];b)). In addition, 
the hole feature clearly splits into two around 9 = 75° 
due to a lifting of the spin degeneracy. 

In a first approach the Fermi surface of graphite has 
been approximated using highly elongated ellipsoids. Ac- 
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cording to the Lifshitz-Onsager relation [53], the funda- 
mental frequencies Bp = hA/2iTe are directly propor- 
tional to the extremal cross sectional areas A of the Fermi 
surface. For the maximal extremal orbits the area is given 
by the intersection of a plane with the ellipse, 



Bfi_ (T) Ai_ 4|| Mx A; 



SWM 



B, 



(X A = nab/^Jsm'^ + [a? /h^) cos^ 



(1) 



where a and b are the semi-major and semi-minor axes 
of the ellipse and 9 is the angle between the magnetic 
field and the c-axis of the graphite crystal. For elongated 
ellipsoids {a/h > 5) this follows very closely nb^/cosO 
for 9 < 60° i.e. follows closely the behavior of a 2D 
cylindrical Fermi surface. 

For our data, matters are further complicated by the 
observed spin splitting for 9 > 70°. In order to include 
spin splitting in our simple model, we note that the os- 
cillatory term can be written as, 



Tosc OC cos 



\Ef\ 



i^T9*f^B)B 




(2) 



where Ef is the Fermi energy, m* is the angle dependent 
effective mass, g* is the Lande g- factor, /is the Bohr 
magnetron and is a phase factor. From Equation [2] 
the frequency in the absence of spin splitting is Bp = 
\Ef\m* /he. Thus, we have a simple relation between 
the frequency (period) of oscillations calculated for the 
ellipsoid and the expected frequencies when spin splitting 
is included. 



1/bI^ = i/BpTgyBmEf\ 



(3) 



so that the expected splitting of the period is simply 
±g* fj,B/'2\Ef\ independent of the angle 9. The Fermi en- 
ergy is E'l^ ~ -0.025 eV for holes and Ef, ^ E'} - 272 ~ 
0.0246 eV for electrons with 72 = -0.0243 eV [3]. A 
reasonable fit to the observed splitting is obtained with 
g* = 2.4 for electrons and g* = 4.0 for holes. The value 
for holes is considerably larger than the value of g* = 2.5 
found for both electrons and holes from magnetotrans- 
port [13] . It is not clear if this is due to the movement of 
the Fermi energy which is not taken into account in our 
analysis, or if g* is really larger for holes. The cross sec- 
tional area of the ellipsoids are obtained by fitting to the 
majority electron and hole frequencies at 9 = 0° and at 
high tilt angles 9 > 70°. The parameters used are sum- 
marized in Table |lj The results of such a fit are plotted 
in Fig(2]as thick solid and dot-dash lines for the major- 
ity hole and electron pockets. The simple model fits the 
experimental data remarkably well, reproducing the ob- 
served angular dependence and the electron and hole spin 
splitting. 

The minority carrier frequencies observed in graphite 
have been reviewed by WooUam [53]. The area of the 
neck orbits can easily be calculated at the iJ-point where 



Maj. hole 
Maj. elec. 

Hole Leg(?) 
Hole neck 



4.7 4.49 40.3 
6.45 6.15 43.1 
~]|2 0.95 17.0 
-0.4J^ 0.41 40.3 



9.0 4.33 

7.0 6.26 
17.8 

98 0.41 



WooIIam[23] 

SWM calculation and WooUam [23] 

TABLE I. Summary of fundamental frequencies and areas 
of the extremal orbits (in units of 10^"^ cm~^) found for the 
ellipsoidal and calculated SWM Fermi surface of graphite. 
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fc^ is linear. 



The Fermi 



the dispersion E = hvf^ 

velocity w/ — \/3eao7o/2/i depends only on 70, whose 
value of 3.15 eV is precisely known from magneto-optical 
data [ini [TTJ |T^. The area of _ff-point neck orbit for 
B ± ab is nkj = irEj/h^vj = 0.43 x lO^^ ^^-^-2 
sponds to a frequency of ~ 0.4 T i.e., the large period 
(1.2 T^^) minority carrier frequency of Ref. [23]. The 
neck orbits have their origin in the two interpenetrating 
hole ellipsoids at the H-point [5]. The size of this or- 
bit is expected to increase rapidly with tilt angle as the 
initially small circular orbit is transformed into a large 
figure of eight orbit encompassing both hole ellipsoids (or 
a single ellipsoid if magnetic breakdown occurs at the H- 
point) [23] . The area of these neck orbits, with and with- 
out magnetic breakdown, have been calculated within our 
simple model using the previously determined parame- 
ters for the majority hole ellipse. The only adjustable 
parameter is the interpenetration of the hole ellipsoids 
which was chosen to have the correct minority carrier fre- 
quency ~ 0.4 T. The calculated period of the neck orbits 
and are shown as dashed lines in Fig|2] The neck or- 
bits have the same frequency the spin split majority hole 
and spin split majority hole harmonic at 9 = 90° and 
so cannot be distinguished. Nevertheless, clear features 
correspond to the neck orbit with magnetic breakdown 
are observed in the data for 75 < 6* < 90°. 

WooUam assigned the minority carrier period of ~ 
1 T^^ to the H-point neck orbits, which in view of our 
results cannot be correct. WooUam's data is plotted as 
symbols in Figj2]ja) and (c) and seems to join up nicely 
with the strong feature at around 50 mT^^ in our data. 
The calculated angular dependence for a neck orbit and 
an ellipsoid are shown as broken lines in Figj2j^c). Clearly, 
the angular dependence corresponds to an ellipsoid rather 
than a neck orbit. The angular dependence for an ellip- 
soid fitted to our data and the data of Woollam is shown 
in Figga). 

Finally, we have calculated the SWM Fermi surface. 
For the diagonalization the SWM matrix is truncated 
to a size of 600 x 600. The magnetic field dependence 
of the density of states at Ef is calculated at fc^ = 
(electrons) and kz = 0.35 (holes) assuming a reasonable 
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FIG. 3. (color online) SWM Fermi surface of graphite along 
the H-K-H edge. 



70 = 3.15 eV 71 = 0.375 eV 72 = -0.0243 eV 
73 = 0.443 eV 74 = 0.07 eV 75 = 0.05 eV 
A = -0.002 Ef = -0.02505 eV 

TABLE II. Summary of the SWM parameters used. 



on graphite reveal the 3D character of the Fermi surface 
of graphite. The Fermi surfaces are closed in all direc- 
tions and well approximated by elongated ellipsoids. Spin 
splitting is clearly observed at high tilt angles and has to 
be included in the analysis in order to extract the correct 
Fermi surface. The SWM parameter 73 is significantly 
larger than previously thought. 

We would like to thank Yu. I. Latyshev for providing 
natural graphite samples. 



Lorentzian broadening of the Landau levels. The Fourier 
transform is then compared with the observed frequen- 
cies for B J- ab. The SWM parameters 70 and 71 are pre- 
cisely known from magnetoptical data[Tni[IIl[I5]. Ej and 
73 are treated as fitting parameters. The hole surface is 
rather insensitive to 73 so that the correct hole frequency 
can be obtained by choosing Ef^ and then the electron 
frequency can be tuned using 73. After a few iterations 
this process converges and the correct electron and hole 
frequencies are obtained with Ef = —0.02505 eV and 
73 = 0.443 eV. The SWM parameters used are summa- 
rized in Table Hll The Fermi surface is then calculated 
by diagonalizing the SWM 4x4 matrix in zero magnetic 
field to calculate the in plane dispersion for = — 0.5 
and looking for the crossing with Ef for angles a = 0— 27r 
in the kx — ky plane. The SWM Fermi surface is shown 
in FigjS] and the calculated cross sectional areas are com- 
pared with the measured dHvA cross sections in Table 
m The good agreement confirms that the diagonalisation 
of the truncated 600 x 600 matrix in magnetic field is 
fully consistent with the results of diagonalizing the 4x4 
SWM matrix in zero field. 

While the calculated Fermi surface consistent with the 
majority electron and hole frequencies it cannot explain 
the observed minority carrier frequency which is well ap- 
proximated by an ellipsoid. Inspecting the SWM Fermi 
surface it can be seen that there is no extremal orbit in 
the vicinity of = 0, so that the frequency should not be 
observed except at high tilt angles, which is indeed the 
case for our data. However, this frequency is very clearly 
seen at 6* = in the data of WooUam. This suggests that 
something is missing from the calculated Fermi surface 
so that a significantly different set of SWM parameters 
might be required. Notably, increasing further the trigo- 
nal warping parameter 73 can generate minority carrier 
pockets. 

In conclusion, angular dependent dHvA measurements 
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